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Approximations of fractional Brownian motion using Poisson processes whose parameter sets 
have the same dimensions as the approximated processes have been studied in the literature. 
In this paper, a special approximation to the one-parameter fractional Brownian motion is 
constructed using a two-parameter Poisson process. The proof involves the tightness and iden- 
tification of finite-dimensional distributions. 
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1. Introduction 

The fractional Brownian motion with Hurst index _ff e (0, 1) is a centered Gaussian 
process = {B^ {t),t > 0} with covariance function 

E[B^(t)B«(s)] = + ~\t~ s\^"). (1.1) 

It follows from (1.1) that B^ is self-similar with index H and has stationary increments. 
Unless H — 1/2 (i.e., B^ is Brownian motion), B^ is not Markovian. Moreover, it is 
known that B^ has long-range dependence if H £ (1/2, 1) and short-range dependence 
ii H E (0, 1/2) (see Samorodnitsky and Taqqu [14]). These properties have made B^ not 
only important theoretically, but also very popular as stochastic models in many areas 
including telecommunications, biology, hydrology and finance. 

Weak convergence to fractional Brownian motion has been studied extensively since the 
works of Davydov [7] and Taqqu [16]. In recent years many new results on approximations 
of fractional Brownian motion have been established. For example, Enriquez [9] showed 
that fractional Brownian motion can be approximated in law by appropriately normalized 
correlated random walks. Meyer, Sellan and Taqqu [13] proved that the law of B^ can be 
approximated by those of a random wavelet series. By extending Stroock [15], Bardina 
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et al. [4] and Delgado and Jolis [8] have established approximations in law to fractional 
Brownian motions by processes constructed using Poisson processes. 

Let {N[t),t> 0} be a standard Poisson process, and for all e > 0, define the processes 
X, = {XS).t>Q} by 



1 f\^i)Nir/e-)^^^ 



£ 

Stroock [15] proved that as e tends to zero, the laws of Xg, converge weakly in the Banach 
space C[0, 1] (i.e., the space of continuous functions on [0,1]) to the law of Brownian 
motion. Delgado and Jolis [8] proved that every Gaussian process of the form 



Xt^ f K{t,s)dBs 
Jo 



where _B is a one-dimensional Brownian motion and K a sufRciently regular deterministic 
kernel, can be weakly approximated by the family of processes 



1 



Their result can be applied to fractional Brownian motion. In addition, Bardina and Jolis 
[2] proved that as £ tends to 0, the family of two-parameter random fields defined by 



Y,is,t)^ rr^^(_i)^(-A.yA)dzd2/ 

JqJo ^ 



(1.2) 



converges in law in the space of continuous functions on [0, 1]^ to the standard Brownian 
sheet. Bardina, Jolis and Tudor showed in [4] that as e tends to 0, the family of two- 
parameter random fields 

t{s,t)^ /Y'l/(s,i,x,y)yi^(-l)^(-/^'^/^)dxd2/ (1.3) 



converges in law to the two-parameter Gaussian process 

W{s,t)= f f f{s,t,x,y)B{dx,dy), 
Jo Jo 

where _B is a standard Brownian sheet and the deterministic kernel 



f{s,t,x,y) = fi{s,x)f2{t,y) (1.4) 

can be separated by the integration variables and satisfies certain conditions. As exam- 
ples, the authors include the fractional Brownian sheet, among others. For more infor- 
mation, see Bardina, Jolis and Rovira [3], where an approximation to the d-parameter 
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Wiener process by a d-parameter Poisson process was provided, and Bardina and Bas- 
compte [1], where two independent Gaussian processes were constructed using a unique 
Poisson process. 

We note that in the serial works [1-4, 8], the dimension of the parameter set is always 
the same for the approximating and the approximated processes. Naturally, we will be 
interested in the problem of whether we can approximate the d-parameter fractional 
Brownian motions by r-paramcter Poisson processes if c^t^ r. The purpose of this paper 
is to study this problem in the case of d = 1 and r = 2. We find that for a special 
deterministic kernel function which cannot be separated with respect to the integration 
variables, the answer is affirmative. Below, we introduce the deterministic kernel function. 

In order to study a non-Gaussian and non-stable process arising as the limit of sums 
of rescaled renewal processes under the condition of intermediate growth. Gaigalas [11], 
page 451, introduced the function h{t,x,y), defined as follows. For x,t>0 and y G M, 

h{t, X, y) = {{t + y)/\0 + x),-{yAO + x)+ 



i: 




-y<x,y< -t, 
-t~y<x<~y,y<~t, 
Mo.x]{u-y)du= -(^ -y, ~y<x,-t<y<0, (1.5) 

< a; < -y, -t<y <0, 
otherwise. 

Note that Kaj and Taqqu [12], page 388, interpreted the function 

Kt{y, x) := {t-y)+Ax- {-y)+ A a; = h{t, x, -x - y) 

in the context of the infinite source Poisson model as a function of the starting time y 
and the duration a; of a session that measures the length of the time interval contained 
in [0, t] during which the session is active. Therefore, h{t, x, y) = Kt{—x — y, x) measures 
the length of the time interval contained in [0, t] during which the session with duration 
X and finishing time —y is active. Obviously, by this definition, h(t,x,y) 7^ if and only 
if ?/ < and —y — x <t, that is, x + y > —t. 
In this paper, we define 

9s (t, x,y) :^ h{t + s, x, y) - h{s, x, y) 

l[o,x]('"-y)dM- / l[o,x](u-2/)dw (1.6) 

-t — s J —s 

— s 



l[o,x](w-2/)du= / l[o^^](w-2/-s)dw = /i(t,a;,y + s) 

for all t > 0, a;>0, y £R and any given s > 0. Then, according to the definition of 
an integral of random measure (see [14], Ghapter 3), we can directly verify that for 

i?e(i,i), 



'Cn / / g,{t,x,y)x''~'W{dx,dy) 



'0 
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JR- 



h{t, x,y + s)x"^'^W{dx, dy) 



(1.7) 



d 




[ h{t,x,y)x"-'^W{dxAy) = B"{t) 



JR_ 



where W{(lx^(ly) is a Gaussian random measure on R_|_ x M with control measure Ax Ay 
(see Section 2 for its definition), Ch = H{2H — 1)(1 — H){3 — 2H) and the notation 

= denotes identification of finite-dimensional distributions. In fact, the last equality is 
taken from Gaigalas [11], page 454, although the constant Ch is omitted in Gaigalas' 
representation. Therefore, 



From the representations (1.2) and (1.7), inspired by (1.3), it seems reasonable that 
the law of the process can be approximated by that of some process similar to 
with kernel gs(t,x,y). 

In this paper we define a sequence of processes {Yn{t),t G [0, l]}„>i as follows: 



for He (1/2,1). Here, {Nn{x, y), {x, y) E M_|- x M_} is a Poisson process with intensity n 
(see Definition 2.1). 

The main purpose of this paper is to show that the law of y„ converges to the law of 
for He (1/2,1). We note that the kernel gs{t,x,y) in (1.9) cannot be separated by 
the arguments {x,y), unlike the kernel function in (1.3). This difference is not trivial. As 
we will see in Remark 3.1, it causes many real technical difficulties. 

The rest of the paper is organized as follows. Section 2 is devoted to introducing the 
necessary definitions, notation and the main result. In Section 3 we prove the family of 
processes {Yn{t)}n>i given by (1.9) is tight in C[0,1]. In Section 4 we prove that as n 
tends to infinity, the finite-dimensional distributions of {F„(t)} converge weakly to those 
of the fractional Brownian motion B^ with HE (1/2,1) and hence {Yn{t)} converges 
weakly in C[0, 1] to the fractional Brownian motion B^ . 

2. Preliminaries 

We now give the definitions of the Brownian sheet and Poisson processes on M x E. Let 
^ be the Borel algebra on R x R. i/ and /x denote a cr-finite measure and the Lebesgue 
measure on R x R, respectively. 

Definition 2.1. Given a positive constant f3 > 0, a random set function N{-) on the 
measure space (R x R, j/) is called the Poisson random measure with density measure 
jiv if it satisfies the following conditions: 




(1.8) 




(1.9) 
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(i) for every Ag ^ with IJ.{A) < oo, N{A) is a Poisson random variable with param- 
eter l3v{A) defined on the same probability space; 

(ii) i/ Ai, . . . , An S ^ are disjoint and all have finite measure, then the random vari- 
ables N{Ai), . . . ,N{An) are independent; 

(iii) ifAi,A2,...€^ are disjoint and A,) <oo, then N A^) = J^^Li N {A^ 
a.s. 



If N{-) is a Poisson random measure with density measure /3fi, then we define 
Nis,t) 



N{[0,s]x [0,t]), s>0,t>0, 

N{[0,s]x[t,0]), s>0,t<0, 

N{[s,0]x[0,t]), s<0,t>0, 

N{[s,0]x [t,0]), s<0, t<0 



and call N = {N{s.,t), {s,t) G M x M} the two-parameter Poisson process with intensity p 
m M X R. 

It is not hard to see that N = {N{s,t), {s,t) e ]R+ x K_} is independent with A^i = 
{A^(s, t), (s, t) € R+ x M+}, which is the ordinary two-parameter process in R+ x M+, and 
for any {s,t) S R+ x M_, {N{s,ty\ has the same finite-dimensional distributions as those 
oi{N{s,\t\)}. 

Definition 2.2. Consider a random set function W{-) on the measure space (R x 
M, ^, v) such that: 

(i) for every Ag with i^iA) < oo, W{A) is a centered Gaussian random variable 
defined on the same probability space with variance 2i/{A); 

(ii) if Ai, . . . , An G ^ are disjoint and have finite measure, then the random variables 
W{Ai), . . . , W{A„) are independent; 

(iii) if Ai,A2,...e^ are disjoint and A,) <oo, thenW{[jZiAi) = J2Zi^i^i) 
a.s. 

We then call W{-) a Gaussian random measure on R X R with control measure v. In 
particular, if W{-) is a Gaussian random measure on R x R with control measure ^fi, 
then we define 



W{[0,s]x[0,t]), s>0,t>0, 

W{[0,s]x [t,0]), s>0, t<0, 

H^([s,0] X [0,i]), s<0, t>0, 

W^([s,0] X [t,0]), s<0, t<0 



Bis,t) 

and call B — {B{s,t), {s,t) e R x R} the two-parameter Brownian sheet m R x 



Similarly, we have that B = {B{s,t), {s,t) S R+ x R_} is independent of Bi = 
{B{s,t), {s,t) G R+ X R+}, which is the ordinary Brownian sheet in R+ x R+, and for 
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any {s,t) G M x R, wc have B{s,t) = B{\s\, \t\). Hence, from (1.7) it is easy to check that 
B"{t) = ^/2C^ / gs{t,x,y)x"-^B{dx,dy). (2.1) 

Let sgn(x) = 1 if x > and sgn(x) = — 1 if x < 0. We have the following conclusion, 
which essentially parallels Bardina and Jolis [2], Theorem 1.1. The proof is omitted. 

Lemma 2.1. Suppose that N = {-/V(s, t), (s, t) G R x M} is a two-parameter Poisson pro- 
cess with intensity 1 in K x M. For any S > and T > 0, let 



y|^(-l)^(^^'^^^'da;dy (2.2) 

for any \u\ < S , \v\ < T. The finite- dimensional distributions of Bn then converge weakly 
to those of a two-parameter Brownian sheet B ~ {_B(u,w), |u| < 5, \v\ < T}. 

Naturally, (2.1) and (2.2) suggest that wc consider the following approximation of B^ 
for i?e (1/2,1): 

no 
g,{t,x,y)x"'^./m\{-l)^"^^'y^ dxdy (2.3) 
-n 

for n g N = {1, 2, . . .} and t e [0, 1], where = {Nn{x, y)} = {N{xy/n, y^/n)} is the two- 
parameter Poisson process with intensity n. 
The main result of this paper is as follows. 

Theorem 2.1. For any fixed s > 0, the law of the process {Yn{t),t e [0, 1]} given by (2.3) 
converges weakly to the law of {B^ {t),t G [0, 1]} for H G (i, 1) in C[0, 1]. 

Remark 2.1. If we define h(t,x,y) :~ h{t,x,^y), then h{t,x,y) has a more natural 
physical interpretation. It measures the length of the time interval contained in [0,t] 
during which the session with duration x and finishing time y is active. Define 



%{t):^n^/2C^ / Ut,x,y)x"-^^{-lf-'^^^y'^dxdy, 
Jo Jo 

where gs{t,x^y) = h(t -\- s,x,y) — h{s,x,y) = h{t,x,y — s). It is then easy to see that 
Theorem 2.1 holds for Y,. 



Note that s is a given positive number. In the sequel we will treat it as a constant. In 
addition, since the parameter Ch cannot affect our discussion, wc will take it to be 1 in 
order to simplify matters. 

We now introduce some auxiliary notation. In the sequel, oo and — cxj will denote 
positive infinity and negative infinity, respectively. For all x' > x > 0, y' > y > 0, by 
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Cairoli and Walsh [6], we define 

A(,,^)iV(x', y') = N{{x, x'] X (y, y']). (2.4) 
For any t e[0,l], x > and y <0, we let 

(j)t{x,y) := gs{t,x,y)/x^~" > 0. 
Define a function Fn,f{x,y) as follows: 



CO ^0 



X Jy 



Fnj{x,y) = ^^\n^ / /(x2,2/2)V^^e-2"W''^-s')-(-^-)?'^ldx2dy2, (2.5) 



where n>0, .i;>0,y<0 and / is a measurable function such that the integral is mean- 
ingful. Obviously, if < / <g, then < Fnjix,y) < Fn,g{x,y). 

3. Tightness of in C[0, 1] 

Let Yn = {Yn{t)^t S [0, 1]} be the process defined by (2.3). The purpose of this section is 
to prove the tightness of the processes {Yn}n>i- 

Proposition 3.1. The family of laws of the processes {Yn}„>i is tight in C[0, 1]. 

To prove the proposition, we need the following lemmas. 

Lemma 3.1. Let Qi = {0 < xi < X2,y2 < 2/i < 0} and ~ {0 < xi < X2,yi < 2/2 < 0}. 
For a non-negative function f{x,y), if J^j'^^ f'^{x,y) dxdy < 00, then for any u,v>0, 



E 



u rO 



/(x,y)v/^(-l)^"(^'^)dxdy 



J -V 

<2(/i(n,/)+/2(n,/)). 



(3.1) 



where 



h{nj)^n' f l[{f{x,,y,)^^\)c-'"^--y---^y^'>dxidyidx2dy2, (3.2) 

f 2 

l2inj)=n^ / n(/(-^"J'')v/^)e"'"'"'^'''^''^"^"'""'^''''i^id2/ida;2dj/2. (3.3) 

Proof. Let /(n, 2:1,2, yi,2) = E[(-l)^"^^''^'^+^^''''^'^]- By Fubini's theorem, the left- 
hand side of (3.1) is equal to 

no fu /"O 
/ / fixi,yi)f{x2,y2)\/xi\yi\V^2\y2\I{n,xi^2,yi,2)dxi dyi da;2d2/2- (3.4) 
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Define fia = {a;i > 2:2 > 0, > t/i > 1/2} and = {xi > X2 > 0,yi < 7/2 < 0}. Then, by 
(2.4), /(n,a;i,2,yi,2) equals 

IE[(-_l')^r.(a;i,|yi|)+Af„(a;2,|y2|)j 

(3.5) 

= E[(— l)'^<°'"'^"^'^i'l^il^+'^"'''')^"^'^^'l^^l^]. 

Note that J2i=i ^io,o)^nixi, \yi\) is equal to the sum of the increments of the Poisson 
process over some disjoint rectangles, and the rectangles which contribute to the value of 
/(n, a;i_2, 2/1,2) are those which appear only once. Since two-parameter Poisson processes 
have independent increments, after some simple calculation, we obtain that on 

I{n,xi.2,yi.2) 

— ]E[(_l)A(o,|„il)Af„(a;i,|y2|)+A(^j j„jl)Ar„(a;2,|y2|)+A(^^,0)Afn(a;2,|yi|)j 

(3.6) 

= E[(-l)'^(o-i«ii)^"(^i'l^"l)]E[(-l)^<-i'i«ii)^"(''^'l^"l^]E[(-l)'^(-i'")^"(''"'l^il)] 
= cxp{-2n(xiyi - X2y2)}- 

Using the same method as above, we obtain that 

-f("-, a;i, 2, 2/1,2) = cxp{-2ri[a;i(y2 - yi) - (a;2 - 3:^1)2/2]} on 172, (3.7) 
-f(",a;i,2,yi,2) = cxp{-2r7,[a;2(2/i - 2/2) - (a;i - .T2)2/i]} on ^3, (3.8) 
^(«,a;i,2,2/i,2) = cxp{-2n(a;22/2 - a;i2/i)} on ri4. (3.9) 

Substituting (3.6)"(3.9) into (3.4) and using a change of the integration variables if 
necessary, we can easily obtain (3.1). □ 

Lemma 3.2. // J^J^^ p{x,y)dxdy <oo, then Ii{n,f) defined by (3.2) is such that 

hinj)<- / f\x,y)dxdy. (3.10) 

o Jo J-00 

Proof. Define B ^ {0 < xi < X2 < 2xi}, C = {2yi < 2/2 < 2/i < 0} and A = S n C. Let 
lf{n), Ii[n) and Ii {n) denote the integral (3.2), where fii is replaced by A, ili\B and 
rii \ C, respectively. Then, 

hin, /) < I^{n, / ) + /f (n, / ) + /f (n, /). (3.11) 

Using the elementary inequality 2ab <a^ + 6^, from (3.2), we have that 

/i^(",/)<i(/ii(n)+/i2(n)), (3.12) 



Approximations of fractional Brownian motion 1203 
where is 

[ /2(a;i,yi).Ti|2;i|e-2"(-iw-x2y2)d^^d2/idx2d2/2 

J A 

no r2xi /.j/i 
/ / /'(xi,yi).Ti|yi|e-2"(-i?'i-=?'^)da;idyidx2dj/2 
-oo J x\ J 2yi 

and Ii2{n) is 

I /'(a;2,y2)a;2|y2|e-'"("^^^-"^^^)dxidyidx2dy2 

J A 

/ / /'(a;2,y2)x2|y2|e-'"("^^^-"^^^Mxidyid.T2d2/2. 

-OO J x\ J 2yi 

Since, for any (xi, yi, X2, 2/2) G ^^i, 

xiVi - X2y2 > {x2 - xi)\yi\ + {\y2\ - |yi|)a;i, (3.13) 

Iii{n) is bounded by 

/ / /'(xi, yi)xi|2/i|e-2"[(-^-i)l?^il+(ly=l-l?'il)-i] dxi dyi da;2 dya- 

Integrating with respect to X2 and then with respect to y2 in the above integral, we 
obtain the fohowing bound: 



-I poo pO 

hi{n)<- / /2(a:i,yi)dxidyi. (3.14) 

Furthermore, in the region A, it foUows from (3.13) that 

xiyi - a;2y2 > {x2 - Xi)|y2 1/2 + (|y2| - |yi|)x2/2. 
By the same argument as above, we then have that 

no 
/2(.T2,y2)d.T2dy2. (3.15) 
-00 

Therefore, from (3.12), (3.14) and (3.15), it foUows that 

r poo pO 

lHn,f)<^ / fix,y)dxdy. (3.16) 

° Jo J~oo 



We now consider I^{n,f). Note that 

2(a;iyi - a;2y2) = 2(a;2 - a;i)|yi | + 2(|y2| - |yi|)a;2 
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> {x2 - xi)\yi\ + (I2/2I - + xi\yi\ 
= {x2 - xi)\yi\ + \y2\x1 

> ^X2\yi\ + \y2\x1 

for any (xi, j/i, X2, 1/2) G f^i \ = {.T2 > 2xi > 0,j/2 < J/i < 0}. Define and /i2(n) as 

/ii(n) and /i2(n), respeetively, with fli\B instead of A. Then, 

-00 J 2x\J —00 

/ / /'(2:2,y2)a;2|y2|e-"(l^^l^^+(^/2)^^I^^IMxidyidx2dy2. 

-00 J 2x\J —00 

Integrating in /ii(n) with respect to X2 and and in Ii2{ti) with respect to Xi and 
respectively, we obtain that 

no 
f{x,y)dxdy. 
-oo 

Since I^{n,f) < ^{Iii{n) +/i2(ri)), we have 



/fK/)<2/ / f'{x,y)dxdy. (3.17) 

Jo J ~oo 

Similarly, for any {xi,yi,X2,y2) G fii \ C = {x2 > a;i > 0, 1/2 < 2j/i < 0}, we have that 
2{xiyi - X2y2) = 2{x2\y2\ - \yi\xi) > 2(\y2\ - \yi\)x2 

> {\y2\ - \yi\)x2 + \y1\x2 > {\y2\ - \yi\)xi + \y1\x2 

> 5l2/2|a;i +a;2|2;l|■ 
Thcrefore, using a similar approach to the one above, we have 

no 
f{x,y)dxdy. (3.18) 
-00 

Combining (3.11) with (3.16)-(3.18), we get (3.10). □ 
From (1.8) and Lemma 3.2, we can immediately get the following corollary. 

Corollary 3.1. For each n>0, h{n,4>t) < j^f^" ■ 

Observe that for every non-negative measurable function f(x,y), by (2.5) and (3.3), 

h{nj)= j dxi f{xi,yi)Fnj{xi,yi)dyi. (3.19) 
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The following lemmas concern l2{n, /). We focus on the case where f = (j)f 

Lemma 3.3. For any < t < 1, let S ~ {x > — s < y < — s} and (f>t{x,y) = 

(f>t{x,y)ls{x,y) . There then exist a non-negative function on {x > 0,y< 0} such 

that for all n> 0, 

F„^^^{xi,yi)lsixi,yi) < ^tixi,yi), (3.20) 
and a positive constant Ki which depends only on H , such that for all n> 0, 



hin, ^t)< dpt{xi,yi)^t{xi,yi) Axi dyi < i^i^ ^-^t'" ■ (3.21) 

Proof. By (1.5) and (1.6), 

tjx^"^ , —y — s < X, y < —t ~ s, 

{x + t + 3 + y)/x^~^, —t — s — y < X < — y — s, y < ^t — s, 

Mx, y) = { -{y + s)/x^-", -~y~s<x, -t-s<y<-s, (3.22) 

1/x^^^ , <x < —y — s, —t — s < y < — s, 

0, otherwise. 

Let Si = S n {0 < X < ~y ~ s} and S2 = Sn{x>-y- s}. Then, 

K,4>t(^^^y^)'^Si{xi,yi) = 'Lsiixi,yi)y/x^{l2i +I22 + I23), (3.23) 

where 

-S-Xi /•-S-J/2 



l2i=n^ / 0t(a;2,y2)v/^7&e-2"[-i(y^-^i)-(-^-i)«^ldx2dy2, 

Jyi Jxi 

l22=n' / 0t(x2,2/2)y^^&e-2"[-ife^-«i)-(-^--i)«^ld.T2d2/2, 



S pOO 



l23=n^ / </'t(a;2,y2)y^^ky^e-2"[-ito^-^i)-(-^-i)^^ldx2dy2. 

For G 5i, from (3.22), wc have that 

l2i=n^ [ ' ' ''^7^/^e-2"[-ito=-yi)-(-^-i>y^Ux2dy2 



^V^"'/ Y '^7^|y2|e-^"[^^(^=-^^)"(^^-^^)''^]dx2dy2 

(3.24) 



<^H-i/2^i±£|^ ' ■'''e-2-ife-^i)dy2 

^ff-1/2 
< C's =:^i(xi,yi), 
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where Cs = ^(1 + s)/s is a constant. Similarly, for {xi,yi) G Si, 



/22 = nM / ^i_^/^e-2"[-^to=-yi)-(-=--i)^=]dx2dy2 

Jyi J-s-y2 X2 

J yi J —s—y2 

(3.25) 



2 



f_S_y,)ff-l/2 

4x1 



and 



/"CXD 



^23="'/ / ^^^^y^e-2"[-ife-^i)-(— i)^^ldx2dy2 

/ ^^|y2|e-'"["^(^^"^^)"("^-"^)^^lda;2dj/2 (3.26) 



1 

In addition, for (xi,j/i) £ 5*2, 



— S /'OO 



/ / ^75^ v/^7^e-2"[-i(«=-^i)-(— 1)^=1 da;2dy2 

(3.27) 

- /2-g =■ V'4(a;i,yi). 

Define \E't(xi,2/i) as 

V^(V'i(a;i,yi) +'02(a;i,yi) +-03(a;i,yi))lsi(a;i,yi) + V'4(a;l,yl)ls2(a;l:yl)■ 
Obviously, ^t{xi,yi) is positive and (3.20) follows from (3.23)-(3.27). Note that 



J s 

< / (f>t{xi,yi)^t{xi,yi)dxidyi. 



With some basic calculations, we obtain the following results: 

(j>t{x,y)y^tpi{xi,yi)dxidyi^ ^ J J 3/2-g '^V^ 4^2 _ i 
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/ (t>t{x,y)./¥['ip2{xi,yi)Axidyi^'^ [ [ ^ '"^ J/2-h dxidyi 

J Si ^ J-t-sJa 

C't^" 



A{2H~1)H' 



0t(x,y)V^V'3(2^i,2/i)dxidyi = — y J ^2-2g Q^^i = g(2ff ^ ' 

0t(a;,y)V'4(a;i,yi)dxidyi = ^ I" |" ^ ^4-2^^ dxi dyi 



-t-sJ-yi-s Xi 



2H 



8{3-2H)H' 

From the above integrals, (3.21) follows with A'l = sh(2H-i) + 8H(3-2g) + 1JP~^- '-' 

Lemma 3.4. For < t < 1 anrf eac/i constant M < —t — s, let G{M) = {—t — s — y < 
x,M <y < —t — s} and (f>t{x,y) = (j)t{x,y)lG{M){x,y) ■ There then exist a non-negative 
function "^tix^y) on {x> 0,y < 0} such that for all n > 0, 

Pn,4>S^i^yi)'^G{M){xi,yi)<'ift{xi,yi), (3.28) 
and a positive constant K2 which only depends on H , such that 



l2in,(pt)< I Mxi,yi)'^tixi,yi)dxidyi<K2\l^-^t^". (3.29) 
Jg{m) V s 

Proof. From (3.22), we have that 

^„,0,(^l-yi)-^21 +/22, (3.30) 

where 

-t — S /'OO ^ 



/2i=n2/ / _^__y^7^e-2"[-ife^-^i)-(-^-i)^^]dx2d2/2, 

122 =n^ r'-'r''^ + + V2 ^^^^-2n;.,fa.-,,)-(..-.0..] dy2. 

For (xi,2/i) e G(M), since |yi| < |Af|, 

/2iVi22< / * dyi r n2|y2|e-2"[-i(y^-^i)-(-^--i)2'^ldx2dy2 

<.;i^i * 
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Therefore, from (3.30) and (3.31), it follows that 



\M\ t 

s Axl^"' 



Let ^t{x,y) ~ y ^-^tx^ ^/2. Then, (3.28) holds. Furthermore, by some basic calcula- 
tions. 



G(M) 



< 



(t>tixi,yi)^tixi,yi) dxi dyi 



< 



\M\ 



\M\ 



(l>t {xi , yi) dxi dyi 

AI J —t— s — yi ZX-y 



s-yi ^-^i 



\M\ 
s 



2H 



4(l-iJ)(3-2iJ) 
+ I dxi 
1 



dxi 



dxi 



t{xi+t + s + yi) 
t{xi+t + s + yi) 



t-s-yi 
-t-s 



t — s—xi 



2x1 



4~2H 



dxi 
dyi 



dyi (3.32) 



' "1 t{xi+t + s + yi) 
2^ 



4(l-i/)(3-2iJ)(2iJ- 1) V s 
Hence, (3.29) holds for K2 = 4(i^h){3-2H){2H-i) 



□ 



Lemma 3.5. For <t <1, M <-2~s, let G{M) = {-t - s - y < x,y < M} and 
(f)t{x,y) = (j)t{x,y)lQ(^j^j'^(x,y). There then exist a non-negative function ^g^tix,y) on 
{x > 0,y <0} such that for all n> 0, 

Fn,^A^i^yi)'^G{M){xi,yi) < ^t{xuyi), (3.33) 
and a constant C > which is independent of t,M and H, such that 



h{n,(t)t)< / <f)tixi,yi)'^tixi,yi)dxidyi < K2Ct^" , 

lG(M) 



(3.34) 



where K2 is the constant in Lemma 3.4- 
Proof. From (3.22), it follows that 



Mx2,y2) n v/^e-2"[-^(^^-^^)-(— dX2 dy2 

-1 .-—I 



'yi ■Jxi 
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3/2-H 



M 



l\yi\ / C-2"[^ite-yi)-(:';2-xi)i;.] j^^^y^ 



(3.35) 



2x1-" JW\ \V2\ 
t 



2x 



2-H 



e-2"=^i|yilna;i|yi| 



d|j/2|e-'""^l^^l 



e^"- dw. 



nxi\M\ ^ 



Let Q{z) = e~^^z J^^ ^^^"^ for z > 1. Then, is continuous on [1, oo) and Q(l) = 0, 
limz_j.oo (3(^) — 1/2. Hence, there is a constant C > 0, which is independent of M,t 
and iJ, such that Q{z) is bounded by C. Note that on G'(Af), because \M\ > 1 and 
xi > —t — s — M > 1, we have nxi |M| > 1 for aU n > 0. Therefore, (3.35) yields 



2x 



2-H 



c. 



Let ^t{x,y) — Ctx /2. By calculations similar to those in (3.32), (3.34) holds for 



7^' _ 1 

"2 4(l-ff)(3-2ff)(2ff-l) 



□ 



Proposition 3.2. For any 0< ^ < 1; there exists a constant K, independent of t hut 
dependent on s, such that 



n / c^t{x,y)^\{-lf-'^^'y^Axdy 



< Kt 



2H 



(3.36) 



Proof. Using the same notation as in Lemmas 3.1-3.5, taking AI = —2 — s and observing 
that from (3.22), (l)t{x,y) = if {x,y) is not in the set 

{{x,y):y<-s,x>0 and x + y > -t - s} ^ S U G{M) U G{M), 

we can rewrite the left-hand side of (3.36) by 



("/7° Mo^,y)V^\i-^f''^^''^dxdy^ 



tix, y) + 0t(a;, y) + 0* (x, y))y^(-l)^"^"'^) Ax Ay 



2-1 



(3.37) 



I{n,(j)t) :=E 

/O J -n 

which is bounded by 

3/(n, 0t) + 3/(n, 0t) + 3/(n, ^t)- (3.38) 
Note that < (pn4'tT4't ^ 4't- Lemma 3.1, Corollary 3.1 and Lemma 3.3 imply that 

/(n,0t)<2/i(n,0t) + 2/2(n,0t) 
^37 ll + s 



< 



2Ki 



(3.39) 



2H 
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and Lemma 3.1, Corollary 3.1 and Lemma 3.4 imply that 



(3.40) 



Furthermore, from Lemma 3.1, Corollary 3.1 and Lemma 3.5, we have 



I{n,cj)t)<2h{n,(f>t) + 2l2{n,4>t)< ( ^ + 2K2C]t 



,2H 



(3.41) 



Therefore, (3.37)^(3.41) yield that 
"111 



I{n,(j)t)<3 



2Ki V(l + s)/s + 2K2 V(2 + s)/s + 2K2C 



,2H 



Taking K = + 2/vi ^(1 + s) /s + 2K2^J{2 + s) /s + 2K2C] then completes the proof 
of Proposition 3.2. □ 

Finally, we prove Proposition 3.1. 

Proof of Proposition 3.1. To prove the tightness of {l^n}n>i in C[0,1], it suffices to 
show that for some ?' > there exist two constants M > and 77 > 1 such that for any 
i,i'e[0,l]. 



g.(t,x,y)-g.(t-,x-,,) ^^_^^^^(,,^)^^^ 



E 

A Jo J-n 

which follows from the criterion given by Billingsley (see [5], Theorem 12.3) and the fact 
that our processes are null at the origin. 

Without loss of generality, let t>t'. Note that from (1.6), we have 

9s{t,x,y)- gs{t',x,y) = h{t + s,x,y) - h{t' + s,x,y) 
= 9t'+s{t-t' ,x,y). 

By Proposition 3.2. it is easy to check that the above inequality holds for r = 2. ij = 
2H>1 and M ^ + 2Ki^/{l + s)/s + 2K2^{2 + s)/s) + 2K2C], which completes 
the proof of Proposition 3.1. □ 

Remark 3.1. Compared with the proofs of tightness in Bardina et al. [2, 4], we have 
found that under the condition that the kernel / can be separated by its arguments 
{x,y), the calculation of I{n,f) is transformed to the calculation of Ii{n,f), which is 
relatively simple; see the proofs of Lemmas 3.1 and 3.2 in [2] and the proof of Lemma 3.1 
in [4]. However, in our case, the kernel / cannot be separated by the arguments {x,y), so 
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we need to discuss /2(n, /). From our proof, we can see that the calculation of l2{n,f) 
is more complicated and delicate than that of Ji(n,/). In addition, the fact that the 
kernel / cannot be separated by the arguments (x, y) also creates some difficulties in the 
identification of the limit law; see the proof of (4.7) in the next section. 

Remark 3.2. By Proposition 10.3 in [10], page 149, Proposition 3.1 also shows that 

P(r„eC[0,i],7ieN) = i. 
4. Limit law of Yn 

In this section, we proceed with the identification of the limit law. We will prove the 
following proposition. 

Proposition 4.1. The finite- dimensional distributions of Yn = {Yn {t),t € [0, 1]} defined 
by (2.3) converge weakly, as n tends to oo, to those of a fractional Brownian motion 
B" = {B"{t),te [0,1]} with Hurst index H e (1/2,1). 

Proof. For each fc e N, ai , . . . , G M and < ti < ^2 < • • • < ifc < 1, we define 



k 



k 




and 




It suffices to prove that for any ^ € R, as n — >■ oo. 



J{n) |E[cxp(iCL„)] - E[exp(i^[/)] | 0. 



(4.1) 



For r > 1 + s, let 




dxdy 



and 




k 



T 



■0 



where B{x,y) is given by Lemma 2.1. Let 



Ji{n,T) 
J2(n,T) 



|E[cxp(iCL„(T))]-E[exp(iet/(r))]|, 

|E[exp(iCL„)]-E[cxp(ie£„(r))]|, 

|E[exp(i^[/(T))]-E[exp(i^f/)]|. 
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Then, 



J{n) < Ji{n,T) + J2{n,T) + MT). 



(4.2) 



Below, we estimate Ji{n,T), J2{n,T) and Jy,{T), respectively. 
(1) We estimate Ji{n,T). 

Noting that (f>tj(x,y) is a non-negative measm'able function on {x>0,y< 0}, we can 
find a sequence of elementary functions g'"'^ (a;, y) such that 

0<q"''Hx,y)<(l)tj{x,y) and q™'^ (x, y) (x, y) a.e. as m -> cx3. (4.3) 

Then, by the dominated convergence theorem, it follows from the fact that J^J^^lfptj {x, 
y)]^ dxdy < cxj that as m — >■ cx), 



T /.O 



/o J-T 
For any m^j,n > 0, define 



{(f>t,{x,y)-q"''\x,y)) dccdy^O. 



(4.4) 



Yn'^^n ff q^^={x,y)^\{~lf-^^^y^dxAy, 

Jo J-T 

Brn.j= r f g"'^(x,2/)B(dx,dy). 

Jo J-T 

By Lemma 2.1, we can readily verify that for fixed m € N, as n — >■ oo. 



Jii{n,T,m) 



E 



cxpl iC^flj^; 



exp ie^a.B™^^ 

V 3 = 1 



•0 (4.5) 



because Y^'^ is essentially a linear combination of increments of i?„ defined by (2.2), 
and B"^'^ is the same linear combination of the corresponding limits of increments of 
We further define 



Ji2(n,r,m) 

Ji3{T,m) 
Then, for any n,TO, 



E[exp(iCL„(T))]-E 
E[exp(i^C/(T))] - E 



expfie^a.r- 

V j=i 

expk^a.S"'- 



Ji{n,T) < Jii{n,T,m) + Ji2{n,T,m) + Ji3{n,T,m). 



(4.6) 
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Below, we will show that for all fixed T, there exists some jri'm) > such that for any 
n > 0, 



Ji2(n, T,m) < 7t(to) 
asm— 5-00. To this end, let j(a;, y) = (/jt^. (x, y) — g'"'^ (x, y). Define 

fni,j{x,y) ■■= fm,j{x,y)llo^T]x[-tj-s,0)ix^y), 
fm,j{x,y) := fr,i,j{x,y)l[o,T]x[-T-t,-s){x,y). 

By (4.3), as m — > oo, 

<^t, ix,y)> fra,j ix,y)^0 a.e. in [0, T] X [-T, 0] . 

Define 



(4.7) 



(4.8) 



Rj{n,m,T) = E 



n / /™,,(x,y)y^(-l)^"(-^'y)dxdy 



J-T 



Then, by Lemma 3.1, 



[Rj{n,m,T)f < n^E 



T M 



(/™.,(2:,2/) + /™,,(x,y))y^(-l)^"(^-^)dxdy 



/O J-T 

<2(/(n, /,„,,) +/(n,/™.j)) 

< 4{Ii{n, frii.j) +hin,frn.j) + hin, fnij) + hinjrnj))- 
Lemma 3.2 and (4.4) show that for any n > 0, as m ^ oo. 



(4.9) 



Il{n,frn,j) < 



37 



/m J (x, y) AxAy=: a-j (m, T) 0. 



(4.10) 



/O J -oo 

Note that fm,j(xi,yi)'^tj{xi,yi) — > a.e. as m— > cx3. From Lemma 3.3 wc know that 
and that 



T rO 



'0 J-tj-s 

< 



f,nj{xi,yi)F^j^^^ {xi,yi)dxi dyi 

'tj{xi,yi)^t, {xi,yi)dxidyi <oo. 



'0 J-tj-s 

By the dominated convergence theorem, as m— > cx). 



T rO 



fni,j{xi,yi)^t, [xi.yi) dxi dyi 0. (4.11) 



J-t,-s 
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In a similar way, we know there are aj{m,T), f3j(m,T) such that for all n > 0, as m — > oo, 
Iiin,fm^j)<aj{m,T)^Q and hinjn,,,) < pj{m,T) ^ 0. (4.12) 
On the other hand, using the mean value theorem, we obtain that 



Ji2(n,T, m) < fc|^| max E 
i<i<fc 



'0 J-T 



= k\^\ max (\a-i\Rj{n,m,T)). 

l<j<k 

Hence, (4.7) follows from (4.9)-(4.13) with 



7t(to) = 2A;|C| max (|ajU/aj(™,T) + ft (m,T) + aj(TO,r) + /3j(m,T)). 



For Ji3(to,T). we apply the mean value theorem again. Then, as m— > oo, 
Ji3(m,T)<P 



< kf max E 

i<j<k 



< k£ max 

i<i<fe 



J-T 



T M 



J-T 
T rO 



hAx,y)-q''''\x,y)) dxdy 



1/2 



•0. 



(4.13) 



(4.14) 



'0 J-T 

From (4.5)-(4.7) and (4.14), we obtain that for any fixed T, as n— >■ oo, 

Ji(n,r)^0. (4.15) 

(2) In a similar way as was used to prove (4.7), there exists some C{T) > such that 
for n> T, as T oo, 



J2(n,T)<C(r)->0. 



(4.16) 



(3) In a similar way as was used to prove (4.14), we obtain that there exists some 
0{T) > such that as T ^ oo, 



JaiT) < 9{T) ^ 0. 



(4.17) 



Therefore, combining (4.2) and (4.15)-(4.17), we can obtain that J{n) — > as n — > oo, 
completing the proof of Proposition 4.1. □ 
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Note that Theorem 2.1 is an immediate result of Propositions 3.1 and 4.1. Therefore, 
the proof of Theorem 2.1 is complete. 

Acknowledgments 

This work was done during our visit to the Department of Statistics and Probability at 
Michigan State University. The authors thank Professor Yimin Xiao for stimulating dis- 
cussions and the department for its good working conditions. This research was supported 
in part by the National Natural Science Foundation of China (No. 10901054). Thanks 
is also due to the anonymous referees for their careful reading and detailed comments 
which improved the quality of the paper. 

References 

[1] Bardina, X. and Bascompte, D. (2010). Weak convergence towards two independent Gaus- 
sian processes form a unique Poisson process. Collect. Math. 61 191-204. MR2666230 

[2] Bardina, X. and Jolis, M. (2000). Weak approximation of tire Brownian sheet from a Poisson 
process in tire plane. Bernoulli 6 653-665. MR1777689 

[3] Bardina, X., Jolis, M. and Rovira, C. (2000). Weak approximation of the Wiener process 
from a Poisson process: The multidimensional parameter set case. Statist. Probab. Lett. 
50 245-255. MR1792303 

[4] Bardina, X., Jolis, M. and Tudor, C.A. (2003). Weak convergence to the fractional Brownian 
sheet and other two-parameter Gaussian processes. Statist. Probab. Lett. 65 317-329. 
MR2039877 

[5] Billingsley, P. (1968). Convergence of Probability Measures. New York: Wiley. MR0233396 
[6] Cairoli, R. and Walsh, J.B. (1975). Stochastic integrals in the plane. Acta Math. 134 111- 
183. MR0420845 

[7] Davydov, Y. (1970). The invariance principle for stationary processes. Teor. Verojatn. Pri- 

men. 15 498-509. MR0283872 
[8] Delgado, R. and Jolis, M. (2000). Weak approximation for a class of Gaussian process. J. 

Appl. Probab. 37 400-407. MR1780999 
[9] Enriquez, N. (2004). A simple construction of the fractional Brownian motion. Stochastic 

Process. Appl. 109 203-223. MR2031768 
[10] Ethier, S. and Kurtz, T. (1986). Markov Processes: Characterization and Convergence. New 

York: Wiley MR0838085 
[11] Gaigalas, R. (2006). A Poisson bridge between fractional Brownian motion and stable Levy 

motion. Stochastic Process. Appl. 116 447-462. MR2199558 
[12] Kaj, I. and Taqqu, M.S. (2008). Convergence to fractional Brownian motion and to the 

Telecom process: The integral representation approach. In In and Out of Equilibrium 

2. Progr. Probab. 60 383-427. Basel: Birkhauser. MR2477392 
[13] Meyer, Y., Sellan, F. and Taqqu, M.S. (1999). Wavelets, generalized white noise and frac- 
tional integration: The synthesis of fractional Brownian motion. J. Fourier Anal. Appl. 

5 465-494. MR1755100 

[14] Samorodnitsky, G. and Taqqu, M.S. (1994). Stable Non-Caussian Random Processes. New 
York: Chapman and Hall. MR1280932 



1216 



Y. Li and H. Dai 



[15] Stroock, D. (1982). Topics in Stochcastic Differential Equations. Bombay: Tata Institute of 

Fundamental Research. MR0685758 
[16] Taqqu, M.S. (1975). Weak convergence to fractional Brownian motion and to the Rosenblatt 

process. Z. Wahrsch. Verw. Gebiete 31 287-302. MR0400329 

Received January 2010 and revised August 2010 



